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Abstract. We propose a discontinuous finite element approximation for a model of 
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1. Introduction 

In this paper we formulate a discontinuous finite element approximation for a model of 
quasi-static growth of brittle fractures in linearly elastic bodies proposed by Francfort and 
Marigo 12]. Their model is based on the classical Griffith's criterion which involves a com- 
petition between bulk and surface energies. To be precise, let ^2 C M'^ be an elastic body, 
be a part of its boundary and let g : be the spatial displacement of Q at 

the points of Od^- According to Griffith's theory, given a preexisting crack Ti C Jl, the new 
crack F and the displacement u : 17 \ F ^ associated to g at the equilibrium minimizes 
the following energy 

(1.1) £{v,g,T):= [ ^uj^vp + |div wj^ + W^(r), 

among all cracks F with Fi C F and all displacements t> : i7 \ F ^ M'^ with v = g on \ F. 
Here Ev denotes the symmetric part of the gradient of v, and TC^ denotes the two dimensional 
Hausdorff measure, while /U and A are the Lame coefficients. The boundary condition is 
required only on \ F because the displacement in a fractured region is assumed not to be 
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transmitted. Let us indicate by S{g, T) the minimum value of among all ?; : \ F ^ 
with V = g on \ T. 

Supposing that the boundary displacement g varies with the time t £ [0, 1], the quasi-static 
evolution t — > T{t) proposed in requires that: 

(1) T(t) is increasing in time, i.e., T(ti) C T{t2) for all < ti < ^2 ^ 1 (irreversibility of 
the process); 

(2) £{g{t),T{t)) < £{g{t),T) for all cracks F such that Us<fr(s) C T (equilibrium condi- 
tion); 

(3) the total energy £{g{t),T{t)) is absolutely continuous in time and (conservation of 
energy) 

^£{g{t),T{t)) = 2fj, [ Eu{t)Eg{t)dx + \ I A\\ u{t) d\\ g{t) dx. 
at Jn Jn 

A precise mathematical formulation of this model has been given by Dal Maso and Toader [S] 
in the case of anti-planar shear in dimension two assuming that the fractures are compact sets 
with a finite number of connected components. Recently Francfort and Larsen ,10 , using the 
framework of 51? y- functions, proved the existence of a quasi-static growth of brittle fractures 
in the case of generalized anti-planar shear and without assumptions on the structure of the 
fractures which are dealt with the set of jumps of the displacements. To be precise, they 
consider as elastic body an infinite cylinder whose section Q C is subject to a displacement 
u S SBV{Q) in the direction orthogonal to fi. The crack at time t on the section is defined 
as 

m := U U {dDnn{u{s) + g[s)))\ , 

s<t 

where Su denotes the set of jumps of u. Moreover the pair (u(t),F(t)) is such that: 

(a) for ah v G SBV(Sl) 

(1.2) I \Vu{t)f dx^n^'^{V(t))< [ \Vv\'^dx + n^-\S^U{dDnn{v^g{t)})ur{t)); 
Jn Jn 

(b) the total energy £{t) := |Vu(t)p dx -\- 7i^^'^{T{t)) is absolutely continuous and 

(1.3) £{t)=£{0) + 2[ [ Vu{T)Vg{T)dxdT. 

Jo Jn 

The aim of this paper is to discretize the model using a suitable finite element method and to 
a give a rigorous proof of its convergence to a quasi-static evolution in the sense of Francfort 
and Larsen. We restrict our analysis to a two dimensional setting considering only a polygonal 
reference configuration C M^. 

The discretization of the domain 0, is carried out, following considering two parameters 
£ > and a g]0, ^[ . We consider a regular triangulation of size e of fl, i.e. we assume that 
there exist two constants ci and C2 so that every triangle T £ contains a ball of radius cie 
and is contained in a ball of radius C2£- In order to treat the boundary data, we assume also 
that is composed of edges of R^. On each edge [x,y] of R^ we consider a point z such 
that z = tx + (l — t)y with t G [a, 1 — a]. These points are called adaptive vertices. Connecting 
together the adaptive vertices, we divide every T E R^ into four triangles. We take the new 



{us'!a '■ i = 0, . . . , Ns} such that uf'a S -4e^a(J7) for all i = 0, . . . , Ng, and such that considering 
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triangulation T obtained after this division as the discretization of 0,. The family of all such 
triangulations is denoted by %^ai^)- 

The discretization of the energy functional is obtained restricting the total energy to the 
family of functions u which are affine on the triangles of some triangulation T{u) £ %_ai^) 
and are allowed to jump across the edges of T(n). We indicate this space by Ae^ai^)- The 
boundary data is assumed to belong to the space AJ-£{0,) of continuous functions which are 
affine on every triangle T € Rg. 

Given the boundary data g G W^''^{[0,1], {n)) with g{t) G AJ'ei^) for all t G [0,1], 
we divide [0, 1] into subintervals of size 6 > fov i = 0, . . . ,Ns, we set gf = g{tf), 

and for all u G Ag^ai^) we indicate by Sj^ (u) the edges of the triangulation T(n) contained 
in on which u ^ gf. Using a variational argument we construct a discrete evolution 

{itffa : i = Q,...,Ns 
the discrete fracture 

r=0 

the following unilateral minimality property holds: 

(1.4) ^ \Vui':f dx< \Vv\^ dx + U Siiv)) \ rf;i-i) . 
Moreover we get suitable estimates for the discrete total energy 

^efa ■= ll^^e,allL2(r2;M2) + (j't'^'a 

The definition of the discrete fracture ensures that T^'^a ^ T^'^a for all i < j, recovering in this 
discrete setting the irreversibilty of the growth given in (1). The minimality property (|1.4)) is 
the reformulation in the finite element space of the equilibrium condition (2). 

In order to perform the asymptotic analysis of the discrete evolution {uf'a : i = 0, . . . , Ng}, 
we make the piecewise constant interpolation in time u^^^^it) = uf'^a and rf^^(t) = rf'^a for all 
tf < t < tf^-^. The main result of the paper is the following theorem. 

Theorem 1.1. Let g G W'^''^{[0,1], H'^ (n)) be such that ||5'(t)||oo < C for all t G [0,1] and let 
ge G W^'^{[0,l],H^{n)) be such that \\geit)\\oo < C, ge{t) G ATe{^) for all t G [0,1] and 

(1.5) ge^g strongly in W^^\[0,l],H\n)). 

Given the discrete evolution {t ai^)} relative to the boundary data gi,, let and £^ ^ 

be the associated fracture and total energy. 

Then there exist 6n — > 0, ^ 0, a„ 0, and a quasi-static evolution {t {u{t),T{t)), t G 
[0,1]} relative to the boundary data g, satisfying (|1.2|) and (|1.3|) . and such that setting Un := 
uil,a^, r„ := Tll^a„, £n ■= £t,a„, the following hold: 

(a) if M is the set of discontinuities ofTl^(r{-)), for all t £ [0, 1] \ we have 

(1.6) Vunit) \/u{t) strongly in L^{n;R'^) 
and 

(1.7) iimn\rn{t)) = n\m); 
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(b) for all t G [0, 1] we have 
(1.8) lim 

n 

We conclude that we have the convergence of the total energy at each time t G [0, 1], and 
the separate convergence of bulk and surface energy for all i G [0,1] except a countable set. 

In order to prove Theorem ll.il we proceed in two steps. Firstly, we fix a and let (5 ^ and 
e -^■ 0. We obtain an evolution {t — > Uaif) : t G [0, 1]} such that Vuf — > Vua{t) strongly 
in L^(r2;M^) for all i up to a countable set and such that the following minimality property 
holds: for all v G SBV{VL) 



(1.9) 



\^ua{t)\'' dx < [ \Vv\''dx + fi{a)nmS,U{dDnn{v^g{t)}))\Ta{t)) , 
Jn 



where /x :]0, ^[^]0, +oo[ is a function independent of e and 6, such that /x > 1, lima^o = 1 
and Ta{t) := Us<t seD '^ua{s) U {do^ n {ua{s) ^ g{s)}). The minimality property ()1.9() takes 
into account possible anisotropics that could be generated as 6 and e ^ 0: in fact, since a is 
fixed, we have that the angles of the triangles in 7^^a(r2) are between fixed values (determined 
by a), and so fractures with certain directions cannot be approximated in length. In the 
second step, we let a — > and determine from {t — > Ua{t) : t G [0, 1]} a quasi-static evolution 
{t — > u{t) : t G [0, 1]} in the sense of Francfort and Larsen. Then, using a diagonal argument, 
we find sequences Sn 0, En ^ 0, and a„ — > satisfying Theorem ll.il 

The main difficulties arise in the first part of our analysis, namely when S,e — > 0. The 
convergence nfa(t) — > Ua{t) in SBV{fl) for i G -D C [0,1] countable and dense is easily 
obtained by means of Ambrosio's Compactness Theorem. The minimality property (|1.9j) 
derives from its discrete version (|1.4j) using a variant of Lemma 1.2 of (lOj : given v G SBV{^}), 
we construct t'f „ G Ae^ai^) such that 



(1.10) ^Ve,a^^v strongly in L^(J^;] 

and 



1.11) hmsup ^1 [(5,.^ U SfJ(*)(<J) \ Tl,{t) 



< 



< f,{a)n^ [{s, u {dDfi n{v^ g{t)})) \ r,(t)] 



where gi{t) := ge{tf) for tf < t < tf^^-^. The main difference with respect to Lemma 1.2 of |1()| 
is that we have to find the approximating functions ^ in the finite element space Ae^ai^)- 
This can be regarded as an interpolation problem, so we try to construct triangulations 
Te G 7^_a(0) adapted to v in order to obtain (|1.10j) and (|1.11|) . In all the geometric operations 
involved, we need to avoid degeneration of the triangles of T(uf ^(t)) which is guaranteed from 
the fact that a is constant: this is the principal reason to keep a fixed in the first step. A 
second difficulty arises when Ua{-) is extended from D to the entire interval [0,1]: indeed it 
is no longer clear whether Vuf ^(t) — > Vua{t) for t ^ D. Since the space Ae^ai^) is not a 
vector space, we cannot provide an estimate on ||Vtif ^(t) — Vu£ ,j(s)|| with s € D and s < t: 
we thus cannot expect to recover the convergence at time t from the convergence at time s. 
We overcome this difficulty observing that Vnf ^ (t) with Ua satisfying a minimality 

property similar to (|1.9jl and then proving Vua = Vua{t) by a uniqueness argument for the 
gradients of the solutions. 
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The plan of the paper is the following. In Section [21 we give the basic definitions and 
prove some auxiliary results. In Section IHl we prove the existence of a discrete evolution. 
In Section 0] we prove the convergence of the discrete evolution to a quasi-static evolution 
of brittle fractures in the sense of Prancfort and Larsen. The proof of minimality property 
()1.9() requires a careful analysis to which is dedicated Section [21 In Section [HI we show that 
the arguments of Section [11 can be used to improve the convergence results for the discrete in 
time approximation considered in jlUj . 

2. Preliminaries 

In this section we state the notation and prove some preliminaries employed in the rest of 
the paper. 

Basic notation. We will employ the following basic notation: 

- is a polygonal open subset of M?; 

- L^(r2;M™') with 1 < p < +00 and m > 1 is the usual Lebesgue space of p-summable 
M™-valued functions, and LP(Q) := LP(0;M); 

- for all /c > 1 and 1 < p < +00, W^'''^{Q) is the usual Sobolev space of functions in 
L^(i7) with distributional derivatives of order 1, . . . ,k in LP(0); we will write H^{Q) 
for W'^^^in); 

- a u G W'''P{0,), Vu is its gradient; 

- is the one-dimensional Hausdorff measure; 

- II • II CO denotes the sup-norm; 

- if / G L2(Q;M"^), II/II denotes the L^-norm of /; 

- for all A C M^, |^| denotes the Lebesgue measure of A\ 

- if /U is a measure on and A is a Borel subset of , ijl\—A denotes the restriction of 
^ to A, i.e. (//L^)(5) := ii{B n A) for all Borel sets S C r2; 

- if fj g]0, -|-cxd[, Oct is such that limCT^o+ e'er = 0. 

Special functions of hounded variation. For the general theory of functions of bounded 
variation, we refer to 4 ; here we recall some basic definitions and theorems we need in the 
sequel. Let A be an open subset of M^, and let u : A ^ M". We say that u G BV{A;W^) 
if u G L^{A]W^), and its distributional derivative is a vector- valued Radon measure on A. 
We say that u £ SBV{A;W^) if n G BV{A;W^) and its distributional derivative can be 
represented as 

Du{A) = / Vu{x)dx+ / {u+{x) -u-{x))(g)iy^dn^-^{x), 
J A J AnSu 

where Vn denotes the approximate gradient of u, Su denotes the set of approximate jumps 
of li, li^ and u~ are the traces of u on Su, i^x is the normal to Su at x, and 7i^^^ is the 
(A'^— l)-dimensional Hausdorfi' measure. The space SBV{A]W^) is called the space of special 
functions of hounded variation. Note that if ti G SBV{A:,W^), then the singular part of Du 
is concentrated on Su which turns out to be countably 'H^~^-rectifiable. 

The space SBV is very useful when dealing with variational problems involving volume 
and surface energies because of the following compactness and lower semicontinuity result 
due to L.Ambrosio (see P, [21, [HI, W)- 



6 A. GIACOMINI AND M. PONSIGLIONE 

Theorem 2.1. Let A be an open and bounded subset o/M^, and let {uk) be a sequence in 
SBV{A;W^). Assume that there exists q > 1 and c > such that 

[ \VUk\Ux+n^-\Su,) + \\Uk\\oo<C 
J A 

for every S N. Then there exists a subsequence (uk^) and a function u S SBV{A;W^) such 
that 

~^ strongly in L^(^;R"); 
(2.1) Vufc^ ^ Vu weakly in L^{A; M^^"); 

In the rest of the paper, we will say that — > u in SBV{A; M") if and u satisfy H2.1() . 
It will also be useful the following fact which can be derived from Ambrosio's Theorem: if 
SBV{A;R'') and if H^-^LSu,, ^ fJ- weakly-star in the sense of measures, then 
n^-^LSu < M as measures. We will set SBV{A) := SBV{A;R). 

Quasi-static evolution of brittle fracture. Let ft be an open bounded subset of with 
Lipschitz boundary, and let dn^l be a subset of dQ open in the relative topology. Let g : 
[0, 1] — > H^{il) be absolutely continuous; we indicate the gradient of g at time t by Vg{t), and 
the time derivative of g at time t by g{t). The main result of (10; is the following theorem. 

Theorem 2.2. There exists a crack r{t) C 17 and a field u{t) S SBV{fl) such that: 

(a) r(t) increases with t; 

(b) u{0) minimizes 

[ iVup dx + n^-^{S^ U{xe dnfl ■ v{x) / c,(0)(x)}) 
Jn 

among all v £ SBV{fl) (inequalities on do^t are intended for the traces of v and g); 

(c) for t > 0, u{t) minimizes 

[ iVvl"^ dx + n^-^ {[S^ U{x£ don v{x) / git){x)}] \ Tit)) 
Jn 

among all v £ SBV{fl); 

(d) Su(t) U {x G : u{t){x) ^ g{t){x)} C r{t), up to a set of Ti'^ ~^ -measure 0. 
Furthermore, the total energy 

£(t) := [ \Vu{t)\^dx + n^-\T{t)) 
Jq 

is absolutely continuous and satisfies 

£{t)=£{0)+2 [ [ Vu{T)Vg{T)dxdT 
Jo Jn 
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for every t G [0,1]. Finally, for any countable, dense set I C [0,1], the crack T{t) and the 
displacement u(t) can be chosen so that 

r(i) = U {S^ir) U {x G don : n(r)(x) / 5(t)(x)}) . 

rG/,T<t 

The main tool in the proof of Theorem 12.21 is the following result [TU', Theorem 2.1], which 
is useful also in our analysis. 

Theorem 2.3. Let ^1 C Jl', with d^l Lipschitz, and let for r = 1, . . . , i (uj^) be a sequence in 
SBV{n') such that 

(1) C n; 

(2) |VuJ^| weakly converges in L^{Q'); and 

(3) lij^ — > strongly in L^{Q'), 

where u' G BV{n') with n^~\Sur) < oo. Then for every <j) G SBV{Q.') with n^-\S^) < oo 
and G -L^(0;M^) for some q G [l,+oo[, there exists a sequence {(j)n) in SBV{^') with 
(pn = <P on Q' \ ^} such that 

(a) (j)n ^ 4> strongly in L^{Q'); 

(b) V(/)„ V0 strongly in L'^{Q.'); and 

(C) W^-l {[S^^ \ Ut=l SurJ \ [S^ \ Ut = l Sur]) - 0. 

/n particular 

(2.2) lim^sup W^-i \ U 5„.^ < W^-i 1^5,^ \ U . 



Hausdorff metric on compact sets. Let A C be open and bounded, and let K,{A) be the 
set of all compact subsets of A. 1C{A) can be endowed by the Hausdorff metric dn defined by 

dniKi, K2) := max < sup dist(x,K2), sup dist(y,-fCi) > , 

[xGKi yeK2 J 

with the conventions dist(x, 0) = diam(j4) and sup0 = 0, so that dni^, K) = if = and 
dj^ (0, K) = diam(^) if / 0. It turns out that K.{A) endowed with the Hausdorff metric is 
a compact space (see e.g. [T^l. 

Triangulations. Let C ]-,g polygonal set and let us fix two positive constants 
< ci < C2. By a regular triangulation of Vt of size e we intend a finite family of (closed) 
triangles Tj such that Q. = Ui^ii Tir\ Tj is either empty or equal to a common edge or to 
a common vertex, and each Tj contains a ball of diameter cie and is contained in a ball of 
diameter C2£- 

We indicate by ^^^(Jl) the family of all regular triangulations of 0, of size e. It turns out 
that there exist < -f?! < 'i?2 < vr such that for all T belonging to a triangulation T G Tl^{Q,), 
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the inner angles of T are between tDi and ■i?2- Moreover, every edge of T has length greater 
than ci£ and lower than C2£- 

Let us fix a triangulation G 7^e(f^) for all e > and let a g]0, |[. Let us consider a new 
triangulation T nested in obtained dividing each T G Rg into four triangles taking over 
every edge [x,y] of T a knot z which satisfies 

z = tx + {1 — t)y, t £ [a,l — a]. 

We will call these new vertices adaptive, the triangles obtained joining these points adaptive 
triangles, and their edges adaptive edges (see Fig.l). 




Fig. 1 



We denote by 7^^a(r2) the set of all triangulations T constructed in this way. Note that for 
all T G Tgr^ai^) there exists < cf < C2 < +00 such that every Tj G T contains a ball of 
diameter cfe and is contained in a ball of diameter C2£- Then there exist < 'df < < tt 
such that for all triangles T belonging to a triangulation T G 7^^a(r2), the inner angles of T 
are between t?f and ??2- Moreover, every edge of T has length greater than cfe and lower than 
cfe. 

We will often use the following interpolation estimate (see (71 Theorem 3.1.5]). If u G 
VF^'^(O) and T G R^, let ut denote the affine interpolation of u on T. We have that there 
exists K depending only on ci,C2 such that 

(2.3) \\UT - u\\iYi,2(T) < K^\W\\w^'^(T)- 

Estimate (|2.3|) holds also for T G Tg^^ai^)- in this case K depends on a. 

Some elementary constructions. The following lemmas will be used in Section 4. 

Lemma 2.4. Let T G %:^a{^), o-nd let I <Z Q be a segment with extremes p, q belonging to 
edges o/T. There exists a polyhedral curve F with extremal points p and q (see Fig. 2) such 
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that r is contained in the union of the edges of those T G T with T n / 7^ 0, and such that the 
following properties hold: 

(1) r = 7j, U 7 U 7ij, where 7 is union of edges of T and 'jp, jq are segments containing p 
and q respectively, and each one is contained in an edge of T; 

(2) there exists a constant c independent of e (but depending on a) such that 

n\T) < cn\l). 




Fig. 2 



Proof. Let {Ti, . . . T^} be the family of triangles in T such that the intersection with / is a 
segment with positive length. For every integer 1 < i < A;, let /, := T,, n / . If /j is an edge of Tj, 
we set Di = Tj. Otherwise let Di be a connected component of Ti \ li such that \Di\ < ^|Tj|. 
We claim that there exists a constant c > independent of e such that 

(2.4) nHdD,) < cn\k). 

We have to analyze two possibilities, namely Di is a triangle, or Di is a trapezoid. Suppose 
that Di is a triangle and that rui is an edge of Di. Let a be the angle of Di opposite to li. It 
is easy to prove that TC^ik) > 7i^{mi) sin a, and so 

n\k) > - sin an\dDi). 
3 

Since < a < 'i?2> sin a is uniformly bounded from below, and hence inequality (|2.4j) follows. 
If Di is a trapezoid, since \Di\ < ^|Tj|, it follows that Tj \ Z?j is a triangle such that its edges 
different from li have length greater than ^cfe. Let a be the inner angle of Tj \ Di opposite 
to li. We have that 

n\li) y^sinacle > i sin a^T^H^A). 
2 2 C2 4 

Since < a < -i?2, inequality follows. 
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By (I^Hl . we deduce that 



moreover, since IJi=i dDi is arcwise connected and contains p, g, we conclude that there exists 
a curve F C {Jl^^ dDi which satisfies the thesis. □ 

Lemma 2.5. There exists a constant c > such that for every segment I C Vl there exists Eq 
with the following property: for every e < Eq, setting TZ{1) := {T G : T Hi ^ 0}, we have 

nHdn{i)) < cn\i). 

Proof Let 7V;(/) := {x £ Q : dist(2;,/) < C2e}. We have that |A/;(/)| = r0-{l)c2e + vrc^e^, and 
hence there exists a positive constant eo such that, for every e < eo, we have that 

|A4(0I < 2n\l)c2e. 

We have that n{l) C Ne{l), and 

4 |A4(0I 



mi) < 



cW 



where jj7^(/) denotes the number of triangles of TZ{1). Then, we have 

n\dn{i)) < 3c2emi) < 3c2£ ^J'^i^^^ < 3cl^2n\i), 

and so the proof is concluded. □ 



A density result. Let A C be open. We say that K <^ A is polygonal (with respect to 
A), if it is the intersection of A with the union of a finite number of closed segments. The 
following density result is proved in fS^. 

Theorem 2.6. Assume that dA is locally Lipschitz, and let u G SBV{A) such that u £ L'^{A), 
Vm G L^(r2;R^), and H^{Su) < +oo. For every e > 0, there exists a function v G SBV{A) 
such that 

(a) is essentially closed, i.e., Tl^{Sy \ Sy) = 0; 

(b) Sy is a polyhedral set; 

(c) V G W''''^{A \ S^) for every k G N; 

(d) \\v - u\\l2(^a) < 

(e) \\Vv - Vu||i2(^.R2) < e; 

(f) \n\Sy)-n\Su)\ <e. 

Let be a relatively open subset of dCl composed of edges lying in d^l. Let us consider 
r^D polygonal open bounded subset of such that $7/) n O = and dQ n d^D = up 
to a finite number of vertices. We set fi' := J7 U 17/) U do^- In Section 4, we will use the 
following result. 
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Proposition 2.7. Given u £ SBV{^1') with u = on ^' \Q and Ti^^'^{Su) < +oo, there 
exists Uh G SBV{^}') such that 

(a) Ufi = in \ 0,; 

(b) Su^ is polyhedral, C Q and Uh G W^'°°{Q,' \ 'S^) for all k; 

(c) Uh ^ u strongly in Lp'iQ.') and Vu^ Vn strongly in L^(r2';M^); 

(d) for all A open subset of Q' with Ti.^{dAr\ Su) = 0, we have 

h ^ 

Proof. Using a partition of unity, we may prove the result in the case :=] — 1, l[x]0, 1[, 
0' :=] — l,l[x] — 1,1[, and do^ ■=] — l,l[x{0}. We set Wh{x,y) := u{x,y — h), and let 
(ffi be a cut off function with if ^ = 1 on ] — 1, l[x] — 1, if^ = on ] — 1, l[x]|, 1[, and 
1 1 V(/?/i| loo < j-^- Let us set Vh '■= (1 — iph)wh- We have that Vh = in ^' \ J7; moreover we have 

Vvh = (1 - (ph)'^Wh - y^hWh- 

Since V(phWh — > strongly in L'^{Q'; M^), we have Vvh —>■ strongly in L?{Vl'; W^). Finally, 
for all A open subset of Q.' with Ti.^{dA n Su) = 0, we have 

limU^iAns^j = n^{AnSu). 

h 

In order to conclude the proof, let us apply Theorem [221 obtaining with polyhedral jumps 
in Q such that Vh G W^^'^{Q.' \ SyJ, \\wh - Vh\\L2{n) + HViu^ - Vvh\\L2{n;m2) < h^ and 
\l-L^~^{S.w^) — 'H^~^{Si,f^)\ < h. If we set := iphg + (1 — ^h)vh: we obtain the thesis. □ 

3. The discontinuous finite element approximation 

In this section we construct a discrete approximation of quasi-static evolution of brittle 
fractures in linearly elastic bodies: the discretization is done both in space and time. 

From now on we suppose that is a polygonal open bounded subset of R^, and that 
C d^l is open in the relative topology. For all e > 0, we fix a triangulation G TZs{i^), 
and suppose that is composed of edges of for all e; we indicate the family of these 
edges by S^. 

We consider the following discontinuous finite element space. We indicate by As^ai^) the 
set of all u such that there exists a triangulation T('u) G T^^ai^) nested in R^ with u affine 
on every T G T(u). For every u G Ae^ai^), we write ||Vn|| for the L^-norm of Vn and we 
indicate by Su the family of edges of T{u) inside across which u is discontinuous. Notice 
that u G SBV{Q) and that the notation is consistent with the usual one employed in the 
theory of functions with bounded variation. Let us also denote by ATei^) the set of affine 
functions in Q with respect to the triangulation R^. Finally, given any g G AJ-£{0,), for all 
u G As^ai^) we set 

(3.1) Sf,{u) ■.= {CeSe : u^g on C}, 

that is Sj^{u) denotes the edges at which the boundary condition is not satisfied. Moreover 
we set 



(3.2) 



S3{u) ■.= SuUS'j,iu) 
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Let now consider g G ^^^'^([O, l];H^{n)) with g(t) G AJ^^{Q) for all t G [0, 1]. Let 5 > and 
let Ns be the largest integer such that S{Ni — 1) < 1; for < i < Ns — 1 we set tf := iS, 
:= 1 and gf := g{tf). The following proposition holds. 

Proposition 3.1. Let e > 0, a €]0, ^] and 6 > be fixed. Then for all i = 0, . . . , Ns there 
exists Ue^a £ •A.e,a{^) such that, setting 



rffa := U 

r=0 

the following hold: 

(a) ||lt£',a||oo ^ ll^f lloo/ 

(b) /or all V £ Ae.ai^) we have 

(3.3) ||V<° f + (5^0(4;°)) < llVt^f + [s^Hv]) 
and 

(3.4) llVnffJP < \\Vvf + ^1 (5^-'(^;) \ rf'^^' 



Proof. The proof is carried out through a variational argument. Let uf'a be a minimum of 
the following problem 

(3.5) min i^WVuf +n\S3o{u))Y 

We set rf'a := 5^o(iif',a). Recursively, supposing to have constructed uf'^a and rf'^a let 
«£ a be a minimum for 



(3.6) min {llVnf + [s^Hu] \ V'^X') } . 

We set rf'a := S^i (ut'^a) U Pffa ^. We claim that problems (|3.5() and ()3.6|) admit a solution 
ui'^a such that Unf'^alloo < 1 1 I loo for all i = 0, . . . , A''^. We prove the claim for problem 1)3. 6() . 



the other case being similar. Let (u„) be a minimizing sequence for problem H3.6|) : since 
is an admissible test function, we deduce that for n large 

Moreover, we may modify ii„ in the following way. If vr denotes the projection in M over the 
interval I := [-jl^f ||oo, \\gi\\ oo], let Un G w4£^a(r2) be defined on each T G T(tt^) as the affine 
interpolation of the values (7r('u„(xi)), vr('u„(x2)), 7r('u„(x3)), where xi, xi and X3 are the 
vertices of T. Note that by construction we have for all n 

ll^tnlloo < ||5,^||oo, l|Vn„|| < ||Vu„||, S^^^{Un) ^ S^^Un), 

SO that (n„) is a minimizing sequence for problem (|3.6jl . We conclude that it is not restrictive 
to assume ||m„||oo < Ibj^lloo- 

Since T(u„) G 7^^a(0), we have that the number of elements of T(?x„) is uniformly bounded. 
Up to a subsequence, we may suppose that there exists an integer k such that T(m„) has 
exactly k elements T^, . . . ,T^. Using a diagonal argument we may suppose that, up to a 
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further subsequence, there exists T = {T^, . . . ,T^} G 7^^a(r2) such that — > T* in the 
Hausdorff metric for all i = 1, . . . ,k. Let us consider T* G T, and let be contained in the 
interior of T*. For n large enough, T* is contained in the interior of and {un)^fi is affine 
with Jrpi I'Vunl'^ dx < C with C independent of n. We deduce that there exists a function 
affine on T* such that up to a subsequence Un ^ u uniformly on T*. Since T* is arbitrary, it 
turns out that u* is actually defined on T* and 

I Vii* P < lim inf / |Vii„Pdx 



Let E Ae^ai^) such that u = on T* for every i = 1, . . . , k: we have 

||Vu||^ < liminf ||Vu„||^. 

n 

On the other hand, it is easy to see that S^^ (u) is contained in the Hausdorff limit of S^^ (un), 
and that 

We conclude that u is a minimum point for the problem H3.6() with ||u||oo < Hfff ||oo- We have 
that point (a) is proved. 

Concerning point (6), by construction we get H3.3() : for i > 1 we have 

llVnffjP + ^1 (5^*(nffJ \rf;r') < ||V«f + (^^f (^) \ rf--'' 

for all V G Ae^ai^), so that 

l|V<if <||V^;f (5^.'(T;)\rf;r'), 

and this proves point (6). □ 

Remark 3.2. For technical reasons due to the asymptotic analysis of the discrete evolution 
ui'^a when 5 ^ 0, e —> and a — > 0, we define uf'^a from Ue,a ^ through problem H3.6() without 
requiring that the adaptive vertices determining rf'_a~^ remain fixed. We just penalize their 
possible changes if they are used to create new fracture: in fact in this case, the surface 
energy increases at each change of a quantity at least of order ae. As a consequence, during 
the step by step minimization, it could happen that some triangles T G 7^^a(r2) contain the 
fracture Pf'^a in their interior. This is in contrast with the interpretation of the triangles 
as elementary blocks for the elasticity problem, but being this situation penalized in the 
minimization process, we expect that it occurs rarely. 

The following estimate is essential for the study of asymptotic behavior of the discrete 
evolution. 

Proposition 3.3. // (uf'^a, rf'^a) fori = 0,...,Ns satisfies condition (b) of Provosition HOI 
setting £e% := UVuf'alP + Ti.^ {^e'^a^ , we have for < j < i < Ns 

i— 1 . „ 

(3.7) < + j^j Vnf';V5(r) dx dr + a' , 
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where 



(3.8) o' 



5 



[4+1 1 fi 

max / \\9{T)\\m{n)dT / \\g{T 
0,...,Ns-lJtS Jq 



r=0,...,Ns-l hs 

Proof. For all < j < — 1, by construction of Ue^f^'^ we have that 

Jf7 



Notice that 



V(5i+i-5')= / '^'v5(T)dr, 



so that 

(3.9) ||V<i+if +7^l(s^l+l(^W)\^f;i) < 



where 

e{5) := max / ||5(T)||/fi(n) 

r=0,...,Ns-l JiS 

From H3.9() . we obtain that for all < j < i < A'^^ 

||v^xffj2 + wi(rf';j < ||vnffj|2 + 7^i(rffj+ 

i-l 



and so the proof of point (c) is complete choosing 

-1 

,<5 



o':=e{6) I \\giT)\\H^n)dT. 







4. The convergence result 



□ 



This section is devoted to the proof of Theorem ll.il As in Sectional let be a polygonal 
open bounded subset of R^, and let C dQ be open in the relative topology. For all e > 0, 
let Re G TZei^) be a regular triangulation of such that is composed of edges of R^. As 
in the previous section, let AJ-£{i^) be the family of continuous piecewise affine functions with 
respect to Rg, and let Ae^ai^) be the family of functions which are affine on the triangles of 
some triangulation T G %^a{^) nested in R^ and can jump across the edges of T. 

In the following, it will be useful to treat points at which the boundary condition is violated 
(see (|3.1() ) as internal jumps. Thus we consider $7/) polygonal open bounded subset of 
such that fi/) n i7 = and D OQd = doVL up to a finite number of points; we set 
fi' := Q U U Given u G Ae^aiS^) and g G ^^^(O), we may extend 5 to a function of 
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H^{^}') and u to a function u £ SBV{i^') setting u = g on ^If). In this way, recalling 1)3.2(1 . 
we have 

S^u) = Su, 

so that the violation of the boundary condition of u can be read in the set of jumps of u. 
Analogously, given u G SBV{Q) and g G H^{Q), we set 

(4.1) S'^iu) := SuU{xe dnn : ^{u){x) + -i{g){x)] 

where 7 denotes the trace operator on d^. We may assume g E H^{Q') using an extension 
operator. We can then consider u £ SBV{0,') such that u = u on 0,, and u = g on Qij. In 
this way we have 

S^{u) = Su up to a set of T^^-measure 0. 
Let us consider g e W'^'^{[0,1], (n)) such that 115(01100 < C for all t G [0,1] and let 
gs e Ty^'H[0, l],H^in)) be such that gelt) G AJ^si^) for all t G [0, 1], 

(4.2) \\9emoo<C 
for all t G [0, 1], and for e ^ 

(4.3) ge^g strongly in W^'\[0,l], H\n)). 

We indicate by {uf'^a, i = 0, . . . ,Ns} the discrete evolution relative to the boundary data g^ 
given by Proposition 13. H and we denote by £^',0. its total energy as in Proposition 13.31 

We assume that g{-) and ge^-) are defined in H^{Q') (we still denote these extensions by g{-) 
and gh{-)), in such a way that ()4.2|) and (|4.3() hold in Q'. Let us moreover set (t) := ge{tf) 
for all tf <t< tf^-^ with i = 0, . . . , A^^ - 1 and (1) := ^^(1). 

Let us make the following piecewise constant interpolation in time: 

4,ait)-=4i f0Vtf<t<tf^, i = 0,...,Ns-l, 

and tif^a(l) := Ue^^ . For all t G [0, 1] we define the discrete fracture at time t as 

s<t 

and the discrete total energy at time t as 

£l,{t) := \\v4^M\' + n' (rf,,(t)). 

We have for all t G [0, 1] 

(4.4) ll<a(i)lloo < Ibf (t)lloc. 

Moreover for all v G Ae^ai^) we have 

(4.5) \\vulMf+n' (5^^'W«,(0))) < iiv^f + ^1 [S^'^('\v)) , 

and for all t g]0, 1] and for all v G ^^^^(il) 

(4.6) l|Vnf,,(t)f < IIVHI' + (5^^'W(z;) \ rf,,(t)) . 
Finally for all < s < t < 1 we have 

(4.7) £l,it) < £l,is) + 2 [ Vul,{T)Vge{T) dx dr + of. 
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where <t < sf < s < and 



(4.8) o; 



max 

r=Q,...,Ns-l JtS 



\9e{T)\\m(n)- 



\9e{r)\\m{n) dr 

For s = we obtain the following estimate from above for the discrete total energy 

(4.9) fj„(t) <4,(0) + 2 r f Vul,{T)VUr)dxdT + ol 

Jo Jn 

where tf <t< tf^^. 

We study the behavior of the evolution {t — > ^(t), t G [0, 1]} varying the parameters in 
the following way. We let firstly e — > and 5 — > obtaining an evolution {t Ua{t), t G [0, 1]} 
relative to the boundary data g with the minimality property (|4.17|) : then we let a ^ 
obtaining a quasi-static evolution of brittle fractures {t u{t) ,t G [0,1]} relative to the 
boundary data g. Finally, by a diagonal argument we deal with {S, e, a) at the same time. 

In order to develop this program, we need some compactness, and so we derive a bound 
for the total energy ^. By (|3.4|) . we have that for all t G [0, 1] 

l|V<,(t)||<||V5f(t)||<(7 

with C independent of 5, e and t. We deduce for all t G [0, 1] 

£ia{t)<£U0) + 2C^ + oi 

Notice that £^ ai^) is uniformly bounded as 5, e vary. Moreover, by (|4.4|) and since ||(?e(t)||oo < 
C for all t G [0, 1], we have that ^(t) is uniformly bounded in L°°{fl) independently of 6,e 
and a. Taking into account (|4.3|) . we conclude that there exists C independent of 5, e, a such 
that for aU te [0,1] 

(4.10) <,(t) + ||<,(t)|U<C". 

Formula ()4.1U|) gives the desired compactness in order to perform the asymptotic analysis of 
the discrete evolution. 

Let now consider 5„ — > and e„ ^ 0: by 1)4. 3() we have 

(4.11) oil ^ 0, 

where of^ is defined in (|4.8j) . By Helly's theorem on monotone functions, we may suppose 
that there exists an increasing function Xa such that (up to a subsequence) for all t G [0, 1] 

(4.12) XnAt) ■■= S^-(^)(nf:,,(s)) j - Xa{t). 
Let us fix D C [0, 1] countable and dense with £ D. 

Lemma 4.1. For all t £ D there exists Ua{t) G SBV{Q) such that up to a subsequence 
independent oft 

^ Mt) m SBV{n). 

Moreover for all t £ D we have 

(4.13) \\Vua{t)f+n'{S^^'Hua{t))) + ||n,(t)|U < C. 
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Proof. Let us consider t ^ D. By (|4.1U|) . we can apply Ambrosio's Compactness Theorem 12. II 
obtaining u G SBV{^) such that, up to a subsequence, uf^ ^(t) — > u in SBV{yi). Let us set 
Ua{t) := n. Using a diagonal argument, we deduce that there exists a subsequence of {5n,£n) 
(which we stiU denote by (5n,en)) such that nf^ ^^{t) Ua{t) in SBV{il,) for all t ^ D. In 
order to obtain inequality H4.13() . we extend uf^ ^(t) and Ua{t) to fi' setting ti^^ ^(t) := gf^{t) 
and tia(i) := fi'(i) on $7/5 ; since Qg^it) — > (7(t) on fi/) strongly in H^{^lj;)), we have that 
^e" a(^) ~^ ^a(i) in SBV{^'), so that we can apply Ambrosio's Theorem, and derive H4.13() 
from (BTTnl) . □ 

The following result is essential for the sequel: its proof is postponed to Section [51 
Proposition 4.2. Let t £ D. For all v G SBV{^1) we have 

(4.14) \\VUa(,t)f <\Nvf+K'^)n\S^'^'\v)\ U S^^'\ua{s))), 

s<t,sGD 

where fi :]0, ^[^]0, +oo[ is such that lima^o /^(o^) = 1- Moreover, Vnf^ ^^{t) \/ua{t) strongly 
in L^{n;R'^). 

We now extend the evolution {t — > Ua{t) : t G D} to the entire interval [0, 1]. Let us set 
for all t G [0, 1] 

s<t,seD 

Lemma 4.3. For every t G [0, 1] there exists Ua{t) G SBV{Q.) such that the following hold: 

(a) for all t G [0, 1] 

(4.15) S^''*\ua{t)) <;LTa{t) up to a set of -measure 0, 
and 

(4.16) \\Vua{t)f + n\S'^^'\ua{t))) + ||n,(t)||oo < C; 

(b) for all V G SBV{n) 

(4.17) \Nna{t)f < \\Vvf + fi{a)n' (^^'^(t;) \ r,(t)) ; 

(c) Vua is left continuous in [0, 1] \ ^ with respect to the strong topology of L'^{Q;M?); 

(d) for all t G [0, 1] \ Ma we have that 

Vuf^^^{t) Vua{t) strongly in L'^{n,R^), 

where Ma is the set of discontinuities of the function Xa defined in H4.12() . 

Proof. Let t £ [0,1] \ D and let tn £ D with tn Z' t. By (|4.13|) . we can apply Ambrosio's 
Theorem to the sequence {ua{tn)) obtaining u G SBV{Q,) such that, up to a subsequence, 
u-aitn) ^ u in SBV{Q). Let us set Ua{t) := u. Let us extend Ua{tn) and Ua{t) to Vl' 
setting Ua{tn) := g{tn) and Ua{t) := g{t) on Q^): we have Ua{tn) — > ?ia(i) in SBV{^'). Since 
W^LS^^^j^) < TL^L-Tait) for all n, as a consequence of Ambrosio's Theorem, we deduce that 
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WiLS'„^(t) < n^LTait). This means L 53^^(0) < W^Lr^Ct), so that KT^ holds. 
Moreover, for all v G SBV{0,), by ()4.14p we may write 

(4.18) ||Vn,(t„)||2 < IIVt;- V5(t)+V5(tn)f + (^'^*^W\ra(tn)) < 

< \\Vv - Vgit) + Vg{tn)f + ^(0)^^ (^S'^^'\v) \ r,(i)) + ^(a)^! (r,(t) \ r,(t„)) , 

so that, since by definition of Ta{t) we have H^{ra{t) \ ra(t„)) — > 0, we obtain that (|4.17j) 
holds; choosing v = Ua{t) and taking the limsup in (|4.18() . we obtain that 

limsup||Vu„(t„)f < \\Vua{t)\\^, 

n 

and so the convergence Vuaitn) Vua{t) is strong in L^(r2,M^). Notice that Vuait) is 
uniquely determined by (|4.15j) and (|4.17j) since the gradient of the solutions of the minimum 
problem 

min |||Vii|p : S^^^\u) C Ta{t) up to a set of "^-measure o| 

is unique by the strict convexity of the functional: we conclude that Vuait) is well defined. 
The same arguments prove that Vn^ is left continuous at all t £ [0,1] \ D. Finally (|4.16j) 
is a direct consequence of (|4.13|) and of Ambrosio's Theorem, and so points (a), (6), (c) are 
proved. 

Let us come to point (d). Let us consider u^^ ^(t) with t Afa, we may suppose that t ^ D, 
since otherwise the result has already been established. By Proposition 14. 21 with D' := D[j{t} 
in place of D, we have that, up to a subsequence, n^^ ^(t) — > u in SBV{Q) such that 

llVnf < \\Vvf + n{a)n^ (^Ss('\v) \ {Ta{t) U ^^^(u))) 

for all V G SBV{n) and Vnf;;_„(t) Vu strongly in L'^{Q.;^^). Let s < t with s G Z); by the 
minimality of ^(s) and by 1)4.10(1 we have 

l|Vnf:,,(s)f < \\Vut^,{t) - Vgtit) + Vg^Mf + KAt) " An,a(5) < 
< l|Vnf:,„(t)f + 2^\\Vgi:{t) - Vgtis)\\ + 

+ ||V5f:(t) - Vgt{s)f + Xn,a{t) - An,a(s). 

Passing to the limit for n — > +00, recalling that ^^^(t) g{T) strongly in H^{Q) for all 
T G [0, 1], we deduce 

\\yua{s)f < llVnf + 2VC^\\Vg{t) - Vg{s)\\ + \\Vg{t) - Vg{s)f + Xa{t) - Xa{s), 

so that, since t is a point of continuity for Aa, Vua is left continuous at t, and g is absolutely 
continuous, we get for s — > t 

||V^x,(t)|p < \\Vu\\\ 
We conclude that Ua{t) is a solution of 

min{||Vf f : S^^^^v) C Ta{t) U S^^^^u) up to a set of -measure 0}, 

so that Vu = Vua{t) by uniqueness of the gradient of the solution. We deduce that 
Vn^^ ^{t) — > Vua{t) strongly in L^(r2;R^), and so the proof is complete. □ 

We can now let a — > 0. 



A DISCONTINUOUS FINITE ELEMENT APPROXIMATION 19 

Lemma 4.4. There exists a„ — > such that, for all t £ D, Ua„{t) u{t) in SBV{i^) for 
some u{t) £ SBV{i}) such that for all v E SBV{^) we have 

(4.19) \\Vu{t)f <\\Vvf +n\S3^'\v)\ J S3^'\u{s))). 

s<t,seD 

Moreover, Vua„{'t) — * Vu{t) strongly in L'^{Q;M?) and 

(4.20) l|Vn(t)f + Wi(5^W(n(t))) + \\u{t)\\^ < C. 

Proof. By (j4.16|) . applying Ambrosio's Theorem to the extensions of Ua{t) to fi' by setting 
Ua{t) := g{t) on Qd, and using a diagonal argument, we find a sequence a„ such that, 
for ah t G D, Ua„{t) u(t) in SBV{n) for some u{t) £ SBV{n) such that K7(H\ holds. 

We now prove that u{t) satisfies property 1)4. 19() . Let v S SBV{Q). Let us fix ti < t2 < 
■ ■ ■ < tk = t with ti G D. We extend v and Ua„{ti) to setting v := g{t) and Ua„{ti) := ^(ti) 
on respectively. Since Ua„(ti) u{ti) in 55^(0') for alH = 1, . . . , A;, by Theorem 12.31 
there exists Vn £ SBV{il,') with Vn = git) on such that Vvn strongly in L^(0';M^) 

and 

(4.21) hm^sup^i (^S.,AUS^.M^ <n' (^^AU^«(*.)) • 
By (|l?THl we obtain 

(4.22) ||Vn,„(Of < WVvnf + Kan)n' (^^-. \ U ' 

so that passing to the limit for n +oo and recalling that fi{a) ^ 1 as a — > 0, we obtain 

\\vu{t)f <\Nvf + n' f 5. \ U • 



1=1 



Thus we get 

\\Vu{t)f < \\Vvf+n^ ^5^W(t;) \ U 5^^*'^(n(ti))j . 

Since ti, . . . ,tk are arbitrary, we obtain (|4.19|) . Choosing v = u{t), taking the limsup in (|4.22jl 
and using K7Tt\ . we obtain Vua„it) Vu{t) strongly in L'^{n;R^). □ 

In order to deal with 6, e and a at the same time, we need the following lemma. 

Lemma 4.5. Let {u{t) : t £ D} be as in Lemma \4-4\ There exist 5n — > 0, e„ — > 0, and 

— > such for all t £ D we have 

nf:,,Jt) ^ n(t) znSBVm . 
Moreover, for all n there exists Bn ^ [0, 1] with < 2^" such that for all t £ [0, 1] \ Bn 

(4.23) ||Vnf:,,Jt)-Vn,„(t)||<i. 
Finally, we have that for all v £ SBV{0,) 

(4.24) l|Vn(0)f + (5f(°)(n(0))) < \\Vvf + (^S3^^\v) 



20 A. GIACOMINI AND M. PONSIGLIONE 

and 

(4.25) StaJO) - \\Vu{0)f+n' {s^^'Hum) ■ 

Proof. Let (a„) be the sequence determined by Lemma 14.41 By Lemma l4.ll for all n there 
exists ((5^, ej^) such that for all t G D and m +oo we have 

and 

V«g,a„(i) ^ Vuajt) Strongly in L\n;R^). 
Moreover by Lemma [4.3l we have that Vu^™ ^ ^"^an quasi-uniformly on [0, 1] as m ^ +oo. 

Let Bn ^ [0,1] with < 2"" such that Vu^n — > Vna„ uniformly on [0,1] \ Bn as 
m ^ +CXD. We now perform the following diagonal argument. Let D = n > 1}. Choose 
nil such that 

\\Vuf^ Jh) - VuaAti)\\ + lln'r Jh) - uMW < 1, 

and 

II ^V' a, (*) - ^^-1 II ^ 1 for ah t G [0, l]\Bi. 
Let m„ be such that 

f^" 1 
llVu,^-,,„(t,) - Vn,„(t,)ll + IWellaM - naAtj)\\ <' for all J = 1, . . . , n 

and 

WVutll^Jt) - VuaAm < \ for all t G [0, 1] \ B,, . 

We may suppose that S%^^ 0, ej^^ 0. Then {6^^, e^^, Un) is the sequence which satisfies 
the thesis. In fact by construction and taking into account (|4.10|) . for all t G D we have 

'^e^",a„ii) u{t) in SBV{Q); moreover the set Bn satisfies (jUSSI)- Notice that uj^" ^^r^ (0) 
satisfies 1)4. 5|l and so ()4.24() and ()4.25|) follow by the F-convergence result of |12j . □ 

Let {5n , £n, dn) be the sequence determined by Lemma 14.51 For all t G [0, 1] let us set 

By Helly's theorem, we may suppose that there exist two increasing functions A and r] such 
that up to a subsequence 

A„ — > A pointwise in [0, 1], 

and 

(4.26) Aa^ T] pointwise in [0,1], 

where Aa„ is defined as in 1)4. 12() . We now extend the evolution {t — > u{t) : t G D} to the 
entire interval [0, 1]. Let us set for all t G [0, 1] 

F(t):= U S^(^)(^(s)), 

s<t,seD 

and let J\f be the set of discontinuities of Tl^{r{-)). Notice that for all t G [0, 1] 

(4.27) 'H\T{t)) < \{t). 
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In fact t G D,\et ti <t2 < . . . < tk = t with G D, consider u;„ G SBV{^}'; M'') defined as 

Wn{x) := {4laAti)ix),...,4l^,^{tk)ix)), 

where we assume that itf^^a„(*j) = 5e"(^«) °^ ^i"- We have ty„ ^ if := {u{ti), . . . ,u{tk)) in 



S BV [Q.' ^M}^) ^ where n(tj) = g{ti) on fi/). Note that for all n we have 5^ = ljf=i -5 t. \ 



so that 

n\s^j < Kit). 

Passing to the limit for n — > +oo and applying Ambrosio's Theorem we get 

(^U = ^'('^-) ^ Iminf WH^^J < A(t); 

we thus have 

^1 |^U59(t.)(^(i,^))^ = nHs^) < x{t) 

and taking the sup over all ti, . . . ,tk, we obtain (|4.27p in D. The case t ^ D follows since 
Tl^{T{-)) is left continuous by definition. 

Lemma 4.6. For every t G [0, 1] there exists u{t) G SBV{Q) such that the following hold: 
(a) for all t G [0, 1] 

(4.28) S'^'W(n(t)) C r(t) up to a set ofTi^ -measure 0, 
and for all t G [0, 1] and for all v G SBV{Q) 

(4.29) l|Vu(t)f < ||V7;f (^^^ (t;) \ r(t)) ; 



(b) Vii is continuous in [0,1] \ (L* UM) with respect to the strong topology of L'^{Q; 



. 



(c) ifj\f is the set of discontinuities of the function rj defined in (|4.26|) . for all t G [0, 1] \J\f 
we have that 

^Ua„{t) Vu{t) strongly in L^(ri,M^). 

Finally 

(4.30) £{t)>£{0) + 2 [ [ Vu{T)Vg{T)dxdT, 

Jo Jn 

where 

Sit) := \\vuit)f+n\rit)). 

Proof. The definition of n(t) is carried out as in Lemma 14.31 considering t G [0, 1] \ D, tn G D 
with tn /" t, and the limit (up to a subsequence) of u(t„) in SBVi^^): 1)4. 28() and (|4.29() 
hold, so that point (a) is proved. It turns out that Vti(t) is uniquely determined and that 
it is left continuous in [0,1] \ D. Let us consider t G [0,1] \ (D U M), and let t„ \ t. 
By Ambrosio's Theorem, we have that there exists u G SBVi^}) with such that, up to a 
subsequence, u(tn) — > u in SBVi^l.). Since i is a continuity point of 'H^(r(-)), we deduce that 
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S^^^\u) C T{t) up to a set of W^-measure 0. Moreover by the minimality property for 
and the fact T{t) C T{t n), we have that for all v G SBV(yi) with 

\\Vu{tn)f < \\Vv -Vg{t) + Vg{tn)f +n' (^^'^(t;) \ r(t„)) < 

< \\Vv - Vg{t) + V5(t„)f + [S^^'\v) \ T{t)) , 

and so we deduce that (|4.29|) holds with u in place of u{t), and that Vn(t„) Vn strongly 
in L^(J7;]R^). We obtain by uniqueness that Vn = \7u{t), and so Vu(-) is continuous in 
[0, 1] \ (Z) U M) and this proves point (6). Point (c) follows in the same way of point (d) of 
Lemma 14. HI 

Let us come to the proof of (|4.,3fl|) . Given t G [0,1] and k > 0, let s'l := for ah 
i = 0, . . . ,k. Let us set li'^(s) := u{s^j^i) for sj^ < s < s^^^^. By (|4.29|) . comparing with 
u{s^_^'^^) — g{s'^_^_i) + g{sf), it is easy to see that 

£{t)>£{0) + 2 [ [ Vu^{T)Vg{T)dTdx + Ok, 
Jo Jn 

where ^ as A; ^ +oo. Since Vn is continuous with respect to the strong topology 
of L^(J7;M^) in [0,1] up to a countable set, passing to the limit for k +oo we deduce 
(IOn|l . □ 

We are now ready to prove the main result of the paper. 

PROOF OF THEOREM im Let D be a countable and dense set in [0, 1] such that e D, 
and let ((5„,e„,a„) and {t u{t) S SBV{Q) : t G [0,1]} be the sequence and the evolution 
determined in Lemma 14.51 and Lemma 14.61 Let us set 

Let AT be the union of the sets of discontinuities of t] and 'H^(r(-)), where r] is defined in 
g^ni)- Let B := flfc^ UhLk^h, where Bh are as in Lemma 1131 since | [JhLk^hl < 2-''+'^, we 
have |i3| = 0. For all t G [0, 1] \ (i3 U Jf) we claim that 

(4.31) Vn„(t)^Vn(t) strongly in ^^(O; M^). 
In fact, since t UhLfc for some k, by Lemma 14.51 we have 

lini||Vnf:,,„(t)-Vn,„(t)|| = 0; 

for t ^ AT, by Lemma 1131 we have that Vua„{t) Vn(t) strongly in L'^{^;R'^) and so 1)4. 31|) 
holds. 

Since gs^ g strongly in ^^-"^'^([0, 1]; H^{Q)), we deduce that for a.e. r G [0, 1] 

V^ejr) ^ Vg{T) strongly in L\n;R^). 

Since £n{0) £{0) by (|4.25|1 and of^ — > 0, by semicontinuity of the energy and by (|4.9j) we 
have that for all t G -D 

(4.32) £{t) < liminf f„(t) < limsup£:„(t) < £{0) + 2 [ [ V u{t)V g{T) dx dr . 

" n Jo Jn 

In view of (|4.30|) . we conclude that for all i G I? 

£{t) =£{{)) + 2 [ [ Vu{T)Vg{T)dxdT, 
Jo Jn 



A DISCONTINUOUS FINITE ELEMENT APPROXIMATION 23 

and since Vn(-) and 7i^{r{-)) are left continuous at t ^ D and so £{■) is, we conclude that 
the equality holds for all t G [0, 1]. As a consequence {t — > u{t) ,t E [0, 1]} is a quasi-static 
evolution of brittle fractures. Let us prove that H4.32() is indeed true for all t G [0, 1]. In fact, 
if t ^ L*, it is sufficient to prove 

(4.33) liminf£:„(t) > 

n 

Considering s > t with s £ D, by (|4.7|) we have 



so that 



liminf£:„(t) > £:(s) - f [ Vu{t)V g{T) dx dr. 
" Jt Jn 



Letting s \ t, since £{■) is continuous, we have (|4.33|) holds. By ()4.32|) we deduce that 
£nit) for all t G [0, 1], so that point (5) is proved. 

We now come to point (a). Since \{t) > H^{T{t)) for all t G [0,1], by (|OT|) and point 
(6), we deduce that A = Tl^(T{-)) in [0, 1] up to a set of measure 0. Since they are increasing 
functions, we conclude that A and Ti^(T{-)) share the same set of continuity points [0, 1] \M, 
and that A = n^{T{-)) on [0,1] \ M. In view of (|4.31|) . point (a) is thus established for all t 
except t G (i3uAA)\AA. In order to treat this case, we use the following argument. Considering 
the measures /i„ := 7i^\—Tn{t), we have that, up to a subsequence, ytn ^J- weakly-star in 
the sense of measures, and as a consequence of Ambrosio's Theorem we have T-0'\—T{t) < 
as measures. Since t ^ M we have /i„(]R^) — > 7i^{r{t)), and so we deduce 'H^Lr(t) = /x. 
Let us consider now Un{t); we have up to a subsequence ii„(t) ^ n in SBV{Q) for some 
u G SBV{n). Setting := c/f;^(t) and u := g{t) on ^d, we have Un{t) — > u in SBV{Q-')^ 

and as a consequence of Ambrosio's Theorem, we get that Ti}\-S3^*\u) < fi = n^\-T{t), that 
is S^^^\u) C r(t). By Theorem 12. 3| we deduce that u is a minimum for 

min{||Vt;||2 : ^^^(t;) C r(t) up to a set of -measure }, 
and by uniqueness of the gradient we get that Vn = Vu{t), so that the proof is concluded. □ 

5. PlECEWISE AfFINE TRANSFER OF JUMP AND PROOF OF PROPOSITION 14.21 

The proof of Proposition 14.21 is based on the following proposition, which is a variant of 
Theorem 12.31 in the context of piecewise affine approximation. 

Proposition 5.1. Given Sn — >■ 0, let g'^ G H^{Q) be such that gl^ G AT£„{^) and g^ — > g'^ 
strongly in H^{Q) for all r = 0, . . . ,i. If Ag^^ai^) is such that —>■ in SBV{^) for 

r = 0,...,i, then for all v G SBV{n) with (^S3\v)j < +oo and Vv G L'^{n;M.^), there 

exists Vn G Ae^^ai^) such that Vn ^ v strongly in L^{0.), Vvn — > Vf strongly in L^(r2;M^) 
and 

(5.1) limsup?^^ (sa'-{vn) \ (j 5^'"«) J < fi{a)n^ Is^^'^v) \ [j 

where fi :]0; ^[-^ with lim„^o+ /^('^) = 1- 

In view of Proposition 15.11 we can now prove Proposition 14.21 
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Proof of Proposition \4.S^ Notice that, in order to prove ()4.14|) . it is sufficient to prove the 
existence of /u :]0; M with hm„^o+ A*!^) = 1 such that, given t G L>, for every = Iq < 
. . . < tr < ■ ■ ■ < ti = t, tr £ D, for all v £ SBV{0,) we have 

(5.2) \\Vuait)f < \\Vvf+^i{a)n' Is^^'Hv) \ (j S^^'^Huaitr))] . 



\ r=0 / 

In fact, taking the sup over all possible to, . . . ,ti, we get 1)4. 14() . 

We apply Proposition O considering gl^ := giliU), g'' := g{tr), < := u^^^^^itr), and 
:= Ua{tr) for r = 0, . . . , There exists /U : ]0; M with lim„^Q+ ^(a) = 1 such that for 

V G SBV{Vl), there exists f„ G Ae^^ai^) with Vvn Vv strongly in L^(r2;M^) and 



limsupW^ S<^'"^^'\vn) \ \J 5^-(*')(txf:,,(t,)) 



< 



r=0 



Comparing n^^ ^ (t) and f „ by means of (|4.6|) , we obtain 

(5.3) \\Vutait)f < \\VVnf+n' (5^-(*)(^n)\rf:,a(t)) < 

so that, passing to the limit for n — > +oo, we obtain that 1)5. 2jl holds. Moreover, we have that 
choosing v = Ua{t), and taking the limsup in (|5.3jl . we get that Vnf^ ^(t) — > Vua{t) strongly 
in L2(J7;]R2). □ 

The rest of the section is devoted to the proof of Proposition 15. 11 It will be convenient, as 
in Section 1^ to consider Qd polygonal open bounded subset of such that fi^) n = and 
90 n dil.D = up to a finite number of vertices; we set := U il/j U do^- We suppose 
that Re can be extended to a regular triangulation of f2' which we still indicate by R^. 

We need several preliminary results. Let us set := — (7^, and let us extend to zero 
on Q/). Similarly, we set := — g^ , and we extend to zero on 

Let 0" > 0, and let C be the set of corners of do^- Let us fix G C M countable and dense: 
we recall that for all r = 0, . . . , i we have up to a set of W^-measure zero 

S,r= (j d*E,,{r)nd*E,,{r), 

where Ec{r) := {x G 11' : z'^{x) > c} and d* denotes the essential boundary (see ^). Let us 
consider 

' 1 
Jj := {x £ U S,r\C : (z') + (x) - (z')-(x) > - for some / = 0, i}, 

r=0 

with j so large that W^dJ^^Q Sz^ \ Jj) < a. Let U he a neighborhood of IJr=o ^^"^ such that 
\U\ < p. Following jlUl Theorem 2.1] (see Fig. 3), we can find a finite disjoint collection of 
closed cubes {Qk}k=i,...,K with center Xk G Jj, edge of length 2rk and oriented as the normal 
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^{xk) to S^r(k) at Xk, such that U^i Qk and Ti^iJj \ [jk=i Qk) < a. Moreover for all 
k = 1, . . . , K there exists r{k) G {0, . . . , z} and ci(r(A;)), C2{r{k)) > such that 



(a) if Xfe G n then Qfe ^ and if x^ G SdO then Qfc H SdO = iJ^, where -ff^ denotes the 
intersection of Qk with the straight line through x^ orthogonal to ^'(xfc); 



(b) n\s^ri,)ndQk) = 0; 

(c) rfe < cH^{S^r(k) n Qjk) for some c > 0; 

(d) < ci(r(fe)) < C2{r{k)) < {z''^^'^)+{x) and C2{r{k)) - ci{r{k)) > ^; 

(e) \ d*E,^^,^k)Mk))] n Qfc) < ark for s = 1,2; 

(f) if s = 1,2, e 5*^;,^(,(fe))(r(A;)) nQfe : dist{y,Hk) > frj) < ar^; 

(g) if := {x e Qk\x ■ v{xk) > 0} and s = 1, 2 

(^■^) l|lE,,(,(fc))(r(fe))nQfc - 1q+ IIlI(Q') < (^^'^fc; 

(h) nH{S^ \ S,r(k)) n Qk) < ark and ^1(5, n dQk) = 0. 



Let us indicate by Rk the intersection of Qk with the strip centered in Hk with width 
2ark, and let us set := {xk ± rke{xk) + sz^(xfc) : s G M} n Rk, where e{xk) is such that 
{e{xk), i^{xk)} is an orthonormal base of with the same orientation of the canonical one. 




and the following hold 




Fig. 3 
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For all B cn', let us set 

TZniB) := {T G : T n B / 0}, T^iB) := {T £ T(z;(^)) : T n B / 0}. 

In the following, we will often indicate with the same symbol a family of triangles and their 
support in R^, being clear from the context in which sense has to be intended. We will 
consider Zn^''^ defined pointwise in \ 5 r{k) and so the upper levels of Zn^'^'^ are intended as 
subsets oi \ S , . 

Lemma 5.2. For all k = 1,...,K there exists G [ci{r{k)),C2{r{k))] such that, setting 
:= {x G TZniQk) ■ Zn^''\x) > c^}, we have 

K 

(5.5) hmsup ^1 ((57e„(Q,)i5^^) \ S K.)) = o., 

fc=l 



and 

(5.6) limsup||l£;fc -Iq+\\l^{q.') < cr'^rl, 

n " 

where 97?,„(Q;.) denotes the boundary operator in TZn{Qk), O'^d Og- — > as a ^ 0. 

Proof. Note that for n large we have Ufc^i^riCQfc) Q U, so that | UfcLi ^n(QA:)| < By 
Holder inequality and since || Vz^H < C for all r = 0, . . . , i, it follows that 

E f |vz;|dx< V||v<||^<(z + i)c'^. 

^=0-^{Ufe7e„{Qfe):r{fc)=r} J J 

Following |in| Theorem 2.1], we can apply coarea- formula for BV- functions (see ^) taking 
into account that Zn^'^^ belongs to SBV{il,') so that the singular part of the derivative is 
carried only by S r{k): since for n large the 7^n(Qfc)'s are disjoint, we obtain 

(5.7) E / ((5^c,n(r(A;)) n 7^„(Qfe)) \ dc < (i + 1)C'^, 
fc=i •^'^ "J 

where Ec^n{r{k)) := {x Q' \ S r{k) : z^n^\x) > c}, and so 

K f'C2(r(k)) I — 

V / 'n'{{dE,^n{r{k))r\nn{Qk)) \ s dc < (i + \)c'^. 

Notice that we can use the topological boundary instead of the reduced boundary of Ec^n{r{k)) 
in 1)5.71) since z^^*^^ is piecewise affine, and so dEc^n{f{k)) \ d* Ec^n{f{k)) / just for a finite 
number of c's. By the Mean Value Theorem we have that there exist G [ci(r(A:)), C2(r(A;))] 
such that 

K 

fc=i 

and taking the limsup for n +oo, we get 1)5. 5() . Let us come to ()5.5() . Since 

^C2(r(fe)),n(?'(^)) ^ £;cfc_„(r(A:)) C £;c^(^(fc))^„(r(/c)), 
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by ()5.4|) we have that for n large 

and so, since \TZniQk) \ Qk\ 0, we conclude that (|5.6|) holds. □ 

Fix k E {1, . . . , K}, and let us consider the family T^[E^). Let us modify this family in 
the following way. Let T G T^{E^)] we keep it if |T n > ^\T\, and we erase it otherwise. 
Let En~^ be this new family of triangles, and let En~ be its complement in T^{lZn{Qk))- 

Lemma 5.3. For all k = 1, . . . , K we have 

K 

(5.8) hmsupj] (dn,,(^Q,)E^n^ \ S^r(u,) = o,, 

K = l 

and 

(5.9) limsup ||lp,fc,+ — ln+l|i ^ 4(T^r|, 

n n ^k 

where Oo- ^ as a ^ ^. 

Proof. Let T £ T^{E^). Since z^*-'^^ is affine on T, it follows that T Ci E^ is either a triangle 
with at least two edges contained in the edges of T or a trapezoid with three edges contained 
in the edges of T. Let 1{T) be the edge inside T where Zn^^^ = c^, where is the value 
determining E!^ (we consider l(T) = if int(T) C E^). In the case T € En~^ as in the case 

T G En , since the angles of the triangles of T{zn^^^) are uniformly bounded away from and 
from vr, arguing as in Lemma [2.4( we deduce that keeping or erasing T, we increase dfi^{^Q^-^E^ 
of a quantity which is less than ctO'iliT)) with c independent of e„. Then we have 

K K K 

Y.nHdn„iQ,)Eii'^\dn^^Q,)E^,)<Y. cWi(/(r)) < ?^H57^„(Q.)^n \ 

k=l k=lT£Tf;(Ek) k=l 

so that taking the limsup for n ^ +oo and in view of 1)5. 5() we deduce that ()5.8|) holds. 

Let us come to (jHUI). Note that \Tj^{dQ+)\ ^ as n ^ +oo. Then if := {T G 

T{zn^^'^) : T C int((5^)}, for n large we have 

IQ+ \ < \At'+ \ E','+\ + \T^HdQt)\ <2\Qt\ E'J + |r„^(5Q+)|, 

where the last inequality follows by construction of En ■ Taking the limsup for n — > +oo, in 
view of ()5.6() we get 

limsup|Q+\^^'+| <2a^rl 

n 

The inequality limsup^ \En~^ \ < 2a^r'^ follows analogously. □ 
For all k = 1, . . . , K and s G M, let us set 

Hk{s) := {x + su{xk), X G Ffc}. 

Lemma 5.4. There exist s^ £]f ''fc) f'^fei (^i^d s~ g] — §?^fc)~f^A;[ such that, setting Hn~^ := 
Hk{Sn) o-'^d Hn~ '■= Hj.{s~) we have for n large enough 

n\H^n^ \ i?^'+) < 20arfc, H^^n'^ n ii;^'+) < 20c7rfc. 
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Proof. By ()5.9|) we can write for n large 



SO that we get G]f ^rfc[ with 

■H\Hk{s+)\E^^^+)<2^ark. 
Similarly we can reason for s~. □ 

Let Sn~^ be the straight line containing Hn^: up to replacing Hn~^ by the connected 
component of Sn~^ r\TZn{Qk) to which it belongs, we may suppose that Hn'"" \ En'"" is a finite 
union of segments with extremes Aj and Bj belonging to the edges of the triangles of 
T^{TZn{Qk)) such that for n large 

By Lemma im for all j there exists a curve inside the edges of the triangles oi'T^{TZn{Qk)) 
joining Aj and Bj and such that 

(5.10) ■H\L+)<cn\l+), 
with c independent of e„. Let us set 

7^'+ := L+ U ^1^2 U L+ U • • • U U L+ . 

Similarly, let us construct 7n' relative to i^n' n E^ . Note that for n large enoug h 7n + n 
Hk{(y) = 0, 7n'~nffA,.(— fj) = 0, and 7n''''n7n'~ = 0. Let us consider the connected component 
of IZniQk) \ 7n'^ containing Hk{a). Similarly, let us consider the connected component 
of TZniQk) \ In containing Hk{—a). For n large enough, by (|5.1()|1 

( ■m— 1 \ m 

i=i J j=i 

A similar estimate holds for d-ji^^^Q^^C^ . 

Let En be the family of triangles obtained adding to En those T G En such that 
T Q , and subtracting those T G En^ such that T Q . Let ~ be the complement of 
Ei'+ in T„^(7^„(Qfc)). 

We claim that there exists C > independent of n such that for all = 1, . . . , X and for 
n large 



(5.12) [dn.i^Q,)K^^ \ dn.(Q,)K'^) < Car,. 

In fact, let ( be an edge of d-ji^(^Qi^^En^\d-ji^(^Q^-jEn~^ , that is ( belongs to a triangle T that has 
been changed in the operation above described. Let us assume for instance that T G En~ and 
T C (:+. If T' is such that T n T' = C, then T' G in fact if by contradiction T' G En'^ , 

then T' G ^n'^ and so we would have ^ ^TZniQk)^^''^ which is absurd. Similarly we get 
T' % . This means that Q ^ (^n„{Qk)^k ^ ^^'^ since the horizontal edges of 7^'^ intersect 
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by construction only elements of En~^ , we deduce that C ^ ^'Rn(Qk)^k \ i^i^i^i-^i)^ ^^'^ W 
we conclude that if^T^ holds. 
We can summarize the previous results as follows. 

Lemma 5.5. For all k = 1, . . . , K there exist two families En^ and En~ of triangles with 
T^iTZniQk)) = Ei'+ U ^n'", Ql\Rk^ E^^ and Q^\Rk<^ E^' , and such that 

K 

(5.13) hmsup^T^i (dnAQ.)^y \ ^.'■(^)) = 

k=l 

where Oo- — > as o" — > 0. Moreover, in the case Xk G d^^l, we can modify En~^ or En^ in 
such a way that En^ (iVt or En ^ ^■ 

Proof. We have that ^^T^ follows from and (f^TT^ . and the fact that Ylk=i^k < c, 

with c independent of a. Let us consider the case Xj. G with \ Rk ^ (the other 

case being similar). From ()5.13p we have that for n large J2f=i {^^'R-niQk)^^ '^ ^ ^fc ) — 

because Zn^'^^ = R£„gh„{r{k)) on and so there are no jumps in . We can thus redefine 
En^ subtracting those triangles that are in obtaining again ()5.13() . □ 

We are now in position to prove Proposition 15.11 

Proof of Provosition \5.li We work in the context of 0,'. For all v £ SBV{il.') with v = 
on H^iSv) < +00 and Vv G L'^{Q';M.'^), we have to construct Vn G SBV{Q') such that 
Vn = Qn on Qd, {vn)\Q G Ae^^ai^), Vn ^ V strongly in L^{0.'), Vvn Vv strongly in 
L2(fi';M2) and 

(5.14) hm^sup^i \ U Sur^ < Ka)n^ (^S, \ (j Sur^ , 

where we suppose that u!^ and are extended to Q' setting nj^ := g^, and u"^ := g^ on 
respectively. 

We set V = g'^ + w, where w G SBV{il') with = on fi/). By density, it is sufficient to 
consider the case w G L°°{Cl'). Up to reducing U, we may assume that ||V5*||/,2(^^.]]j2j < c 
and ||Vt(;||i2((7.]g2) < a. Let i?'^ be a rectangle centered in Xk, oriented as Rk, and such that 
R'j^ C int-Rfc and H^iSw n {Rk \ i?^)) < OTfc. We claim that there exists w^- G SBV{^1') with 
Wa- = w on R'j^ and = in il/j such that 

(1) \\w — Wf^W + \\Vw — VWf^W < a; 

(2) n\S^^ n {Qk \ R'k)) < o^rk, with ^ as a ^ 0; 

(3) nHs^^ \ uf=i Rk) < n^s^ \ uf=i Rk) + ^; 

(4) Sw„ \ UfcLi ^k is union of disjoint segments with closure contained in \ UfcLi ^k] 
2,oo n \ /'l 1^ 



(5) Wa is of class W^'°° onn\[ \Jf^^ Rk U S, 
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In fact, by Proposition 12.71 there exists Wm £ SBV{W) with Wm = in 17' \ such that 
Wm — >■ w strongly in L'^{Q'), Vwm — > strongly in L'^{Q';'M?), is polyhedral with 

S;Z^^,w^ is of class W^'^ on 17 \ (uf=i Rk U ^) , and lim^ n^AnS^J = n\AnS^) 
for all A open subset of Q' with '^^((JA n 5^^,) = 0. It is not restrictive to assume that 
l-0-{Sw n dRk) = and TO-{Sw^ n dRk) = for all m. Let V'fc be a smooth function such 
that < Vfc ^ 1) V'fc = 1 on -Rfc and Vfc = outside ii^. Setting V := Z^fcLiV'fc) let us 
consider Wm '■= i'w + (1 — ip)'Wm- Note that Wm — > strongly in VtLim — ^ 

strongly in L^(J7';M'^), = in Qd- Moreover, by capacity arguments, we may assume 
that \(Jfc=i Rk is a finite union of disjoint segments with closure contained in 17\Ua;=i ^k- 
Finally, for m — > +oo, we have 

K K 

n\s^^ \\jRk)^ n\s^ \ U Rk)), 

k=l k=l 



K K 

n\s^^ n U (Qfc \ Rk)) ^ n\Su, n |J (q^ \ Ru)) 

k=l k=l 

and ImisvL'p^l-O' {Sw^^^r\{Rk\R'j.)) < 21-0'{Sw^{Rk\R'k)) ^ ^orfc. Then we can take := Wm 
for m large enough. 

Let Syj^ \ Qk '■= Ujli where, by capacity arguments, we can always assume that 
Ij are disjoint segments with closure contained in \ [jj^^i Qk- We define a triangulation 
Tn £ %n,ai^') specifying its adaptive vertices as follows. Let us consider the families TZn{Qk) 
and Tln{lj) for k = 1,. . . ,K and j = 1, . . . ,m. Note that for n large enough, TZniQki) H 
TlniQk,) = for ki ^ k2, Uniln) n Unilj^) = for ji / j2, and 7^n(Qfc) n Unilj) = 
for every k,j. We consider inside the triangles of TZniQk) the adaptive vertices of T(z^^'^^). 
Passing to TZn{lj), by density arguments it is not restrictive to assume that Ij does not pass 
through the vertices of Re„ and that its extremes belong to the edges of R^^. Let ( := [x,y] 
be an edge of TZn{lj) such that Ij f\ C, = {P}- Proceeding as in ^21; we take as adaptive 
vertex of C, the projection of P on {tx + (1 — t)y : t G [a, (1 — a)]}. Connecting these adaptive 
vertices, we obtain an interpolating polyhedral curve Ij with 

(5.15) li}{lj) < fi{a)n\lj), 

where /i is an increasing function such that lima^o /^(o) = 1- Finally, in the remaining edges, 
we can consider any admissible adaptive vertex, for example the middle point. 

Let us define Wn G SBV{^1') in the following way. For all Q^, let Wn be equal to Wa on 
T^n{Qk)\Rk, equal to the reflection of Wa^Q+\^j^^ with respect to Hk{a) on En'^nRk and equal 

to the reflection ofwa^Q-x^^^ with respect to Hk{—a) on En~nRk, where En^ are defined as in 
Lemma 15.51 On the other elements of T„, let us set Wn = w^- Notice that Wn = on and 
that inside each TZniQk)-, all the discontinuities of Wn are contained in d-ji^i^Q^^-^En^ UVkU P^^, 
where P^^ is the union of the polyhedral jumps of Wa- in TZniQk) and of their reflected version 
with respect to Hk{±a). By Lemma 1^31 and since ^^^^^'^^(Vfc U P^^) < Oq- with Oo- —> as 
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o" — > 0, and (^Ur=o \ UfcLi Qkj < 2c7, we have that 



r=0 



< \^S^^ \ U Qfcj + lim^sup \ U %J n nniQk)j < 



r=0 / V=0 k=l 



\ r=0 / 

and since || Vzi;CT||j;^2((7.iR2) < we get for n large 

We now want to define an interpolation Wn of Wn on T^. Firstly, we set Wn = on all regular 
triangles of VId- Passing to the triangles in lZn{Qk) (see fig. 4), by Lemma 12.51 we know that 
for n large enough, we have 

n\dnn{v,)) < cn\Vk), n\dnn{PtJ) < cn\Pl), 

with c independent of n. If T G 7^„(Vfc) U7^„(P^_^), we set Wn = on T; otherwise, we define 
Wn on T as the affine interpolation of Wn- 

Since Vwn is uniformly bounded on TZn{Hk{±cr)), \TZn{Hk{±cr))\ — > and since Wn is 
uniformly bounded in W'^'°° on the triangles contained in lZn{Qk) \ T^n{Vk U U Hk{±cr)) 
we have by the interpolation estimate 1)2.3^ and by H5.16|) 

(5.17) lim^sup||Vu;„f^,^^, ^^^^^^^^ < o.. 
Moreover we have 

K / k \ 

(5.18) limsup^H^ (S^^ \ U 5,.) n7^„(Qfe) < o,. 

k=l \ r=0 / 

Let us come to the triangles not belonging to TZniQk) ior k = 1, . . . , K. For all j = 1, . . . , m, 
we denote by TZn{lj) the family of regular triangles that have edges in common with triangles 
oiTZn{lj)- For n large we have that 'R,n{lji)r\TZn{lj2) = for ji 7^ j2- On every regular triangle 
T Ufc=i T^niQk) U UJLi T^nilj), we define as the affine interpolation of w^- Since Wo- is 
of class Ty2,oo 2^ J. 

is regular, we obtain by the interpolation estimate (|2.3|) 

(5.19) llzDn - ■"^(7 1 1 ^1,2(7.) < Ken\\wcr\\w^-°°- 

Let us consider now those triangles that are contained in the elements of UjLi^nl^j)- Fol- 
lowing [12^, we can define Wn on every T in such a way that Wn admits discontinuities only 
on Ij, and ||Vt(;„||/^cx)(2-) < HVwo-Hoo- Since \TZn{lj)\ as n ^ 00, we deduce that 

(5.20) lim||Vu;„||^,(^^(^^,)^ = 0. 




We are now ready to conclude. Let us consider Wn S Ae^^ai^) defined as Wn ■= On + Wn- 
We have Wn ^ + strongly in L'^{n'). By (EHZI), (j^TT^ . (|^^ we get 

limsup ||Vu)„||^ < llV^i* + Vwo-ll^ + Oct, 

n 

while by ((??TH|) and (|^^^ we have 

limsupH^ (s^„\(j < f,{a)n^ [ \ U sA + o^. 

\ r=0 J V r=0 / 

Letting now a — > 0, using a diagonal argument, we conclude that Proposition 15. II holds. □ 



6. Revisiting the approximation by Francfort and Larsen 

Li this section we show how the arguments of Section|l]may be used to deal with the discrete 
in time approximation of quasi-static growth of brittle fractures proposed by Francfort and 
Larsen in ^Uj . More precisely, we prove that there is strong convergence of the gradient of the 
displacement (in particular convergence of the bulk energy) and convergence of the surface 
energy at all times of continuity of the length of the crack; moreover there is convergence of 
the total energy at any time. 
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We briefly recall the notation employed in ^U]. Let /qo be countable and dense in [0,1], 
and let := {0 = tg < . . . < = 1} such that (/„) is an increasing sequence of sets whose 
union is loo- Let Q Q be a Lipschitz bounded domain, and let dQ = U dQj, where 
SOj is open in the relative topology. Let Q' C 'K^ be open and such that Q, C 0', and let 
g G VF^'^([0, 1]; {Q')). At any time t^, Francfort and Larsen consider u'^ minimizer of 



\J Su^^udnf 

o<i<fc-i 



in {v G SBV{n') : v = g{tl) in \ O}. Setting n"(t) := for t € and r"(t) := 

[Js<t,sein -^f^Cs) ^ ^^/' P'^°^*^ 

(6.1) £''{t)<£^{0) + 2 [ W{T)Vg{T)dxdT + On, te[t^,tl^,[, 

Jo Jq 

where £:"(t) := |Vu"(t)p + W^-^ (r"(t)) and o„ ^ as n ^ +oo. Using Theorem lO 
they obtain a subsequence of (u"(-)), still denoted by the same symbol, such that u^{t) u{t) 
in SBV{9.') and Vn"(t) ^ Vn(t) strongly in L'^{Q.';R^) for ah t G /oo, with u(t) a minimizer 
of 

/ \vv\^dx + n^~^{s.,\v{t)), 

where r(i) := IJse/oo s<t^u{s) U (9il/. The evolution {t ^ u{t), t £ Iqo} is extended to the 
whole [0, 1] using the approximation from the left in time. 

We can now use the arguments of Section |11 Following Lemma 14.61 it turns out that for 
all t G [0, 1] 

(6.2) £{t)>£{0)+2 [ [ Vu{T)Vg{T)dxdT. 

Jo Jq 

Moreover, by the Transfer of Jump and the uniqueness argument of Lemma ESI we have that 
Vn"(t) — > Vii(t) strongly in L^(ri';M^) for all t M, where M is the (at most countable) 
set of discontinuities of the pointwise limit A of 7i^~^{T{-)) (which exists up to a further 
subsequence by Helly's Theorem). Then we pass to the limit in 1)6. 1|) obtaining 

£{t) <£{<d) + 2 ( [ Vu{T)Vg{T)dxdT; 
Jo Jn 

moreover, following the proof of Theorem ll.il we have that for all t G [0, 1] 

£{t) < liminf£:„(t) < limsup.?„(t) = ^(0) + 2 [ [ Vu{t)V g{T) dx dr, 
" n Jo Jn 

and taking into account (|6.2j) we get the convergence of the total energy at any time. Since 

Vn"(t) Vu{t) strongly in L^{n';R^) for every t G loo, we deduce that A = n^-^{T{-)) on 

loo, so that the convergence of the surface energy holds in loo- The extension to the continuity 

times for 7i'^~^(r{-)) follows like in the final part of the proof of Theorem ll.il 
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